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Quantum number-path entanglement is a resource for super-sensitive quantum metrology and in 
particular provides for sub-shotnoise or even Heisenberg-limited sensitivity. However, such number- 
path entanglement has thought to have been resource intensive to create in the first place — typi¬ 
cally requiring either very strong nonlinearities, or nondeterministic preparation schemes with feed¬ 
forward, which are difficult to implement. Very recently, arising from the study of quantum random 
walks with multi-photon walkers, as well as the study of the computational complexity of passive 
linear optical interferometers fed with single-photon inputs, it has been shown that such passive 
linear optical devices generate a superexponentially large amount of number-path entanglement. A 
logical question to ask is whether this entanglement may be exploited for quantum metrology. We 
answer that question here in the affirmative by showing that a simple, passive, linear-optical in¬ 
terferometer — fed with only uncorrelated, single-photon inputs, coupled with simple, single-mode, 
disjoint photodetection — is capable of significantly beating the shotnoise limit. Our result implies 
a pathway forward to practical quantum metrology with readily available technology. 


Ever since the early work of Yurke & Yuen it has been 
understood that quantum number-path entanglement is 
a resource for super-sensitive quantum metrology, allow¬ 
ing for sensors that beat the shotnoise limit [1, 2]. Such 
devices would then have applications to super-sensitive 
gyroscopy [3] , gravimetry [4] , optical coherence tomogra¬ 
phy [5], ellipsometry [6], magnetometry [7], protein con¬ 
centration measurements [8], and microscopy [9, 10]. This 
line of work culminated in the analysis of the bosonic 
NOON state ((|A^, 0) + |0, Y))/v/2, where N is the total 
number of photons), which was shown to be optimal for 
local phase estimation with a fixed, finite number of pho¬ 
tons, and in fact allows one to hit the Heisenberg limit 
and the Quantum Cramer-Rao Bound [11-14]. 

Let us consider the NOON state as an example, where 
for this state in a two-mode interferometer we have the 
condition of all N particles in the first mode (and none 
in the second mode) superimposed with all N particles 
in the second mode (and none in the first mode). While 
such a state is known to be optimal for sensing, its gener¬ 
ation is also known to be highly problematic and resource 
intensive. There are two routes to preparing high-NOON 
states: the first is to deploy very strong optical nonlin¬ 
earities [15, 16], and the second is to prepare them us¬ 
ing measurement and feed-forward [17-19]. In many ways 
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then NOON-state generators have had much in common 
with all-optical quantum computers and therefore are 
just as difficult to build [20]. In addition to the com¬ 
plicated state preparation, typically a complicated mea¬ 
surement scheme, such as parity measurement at each 
output port, also had to be deployed [21]. 

Recently two independent lines of research, the study 
of quantum random walks with multi-photon walkers in 
passive linear-optical interferometers [22-24], as well as 
the quantum complexity analysis of the mathematical 
sampling problem using such devices [25, 26], has led to a 
somewhat startling yet inescapable conclusion — passive, 
multi-mode, linear-optical interferometers, fed with only 
uncorrelated single photon inputs in each mode (Eig. 1), 
produce quantum mechanical states of the photon field 
with path-number entanglement that grows superexpo¬ 
nentially fast in the two resources of mode and photon- 
number [27]. What is remarkable is that this large de¬ 
gree of number-path entanglement is not generated by 
strong optical nonlinearities, nor by complicated mea¬ 
surement and feed-forward schemes, but by the natural 
evolution of the single photons in the passive linear op¬ 
tical device. Whilst such devices are often described to 
have ‘non-interacting’ photons in them, there is a type 
of photon-photon interaction generated by the demand 
of bosonic state symmetrization, which gives rise to the 
superexponentially large number-path entanglement via 
multiple applications of the Hong-Ou-Mandel effect [24] . 
It is known that linear optical evolution of single pho¬ 
tons, followed by projective measurements, can give rise 
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to ‘effective’ strong optical nonlinearities, and we conjec¬ 
ture that there is indeed a hidden Kerr-like nonlinearity 
at work also in these interferometers [28]. Like boson¬ 
sampling [25], and unlike universal quantum computing 
schemes such as that by Knill, Laflamme, and Milburn 
[29], this protocol is deterministic and does not require 
any ancillary photons. 

The advantage of such a setup for quantum metrol¬ 
ogy is that resources for generating and detecting single 
photons have become quite standardized and relatively 
straightforward to implement in the lab [30-36]. The 
community then is moving towards single photons, lin¬ 
ear interferometers, and single-photon detectors all on a 
single, integrated, photonic chip, which then facilitates a 
roadmap for scaling up devices to large numbers of modes 
and photons. If all of this work could be put to use for 
quantum metrology, then a road to scalable metrology 
with number states would be at hand. 

It then becomes a natural question to ask — since 
number-path entanglement is known to be a resource 
for quantum metrology — can a passive, multi-mode 
interferometer, fed only with easy-to-generate uncorre¬ 
lated single photons in each mode, followed by uncor¬ 
related single-photon measurements at each output, be 
constructed to exploit this number-path entanglement for 
super-sensitive (sub-shotnoise) operation? The answer is 
indeed yes, as we shall now show. 

The phase-sensitivity, Aip^ of a metrology device can 
be defined in terms of the standard error propagation 
formula as. 


Aip = 




d{0) 

dip 


( 1 ) 


where (O) is the expectation of the observable being mea¬ 
sured and (f is the unknown phase we seek to estimate. 

The photons evolve through a unitary network accord¬ 
ing to UalU^ = Uija}j. In our protocol, we construct 
the n-mode interferometer tj to be, 

IJ (2) 


which we call the quantum fourier transform interferom¬ 
eter (QuFTI) because V is the n-mode quantum Fourier 
transform matrix, with matrix elements given by. 


M _ _ 


—2ijk7T 
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( 3 ) 


and 0 are both diagonal matrices with linearly increas¬ 
ing phases along the diagonal represented by. 


^j,k = exp i{j - l)(p 
©i.fe =exp 


( 4 ) 


where (f is the unknown phase one would like to measure 
and 0 is the control phase. 0 is introduced as a reference. 


which can calibrate the device by tuning 0 appropriately. 
To see this tuning we combine and 0 into a single 
diagonal matrix with a gradient given by. 


^j,k ’ ^j,k 


Sj^k exp 




( 5 ) 


The control phase 0 can shift this gradient to the optimal 
measurement regime, which can be found by minimizing 
Aip with respect to n and ip. Since this is a shift according 
to a known phase, we can for simplicity assume (and 
without loss of generality) that ip is in the optimal regime 
for measurements and ^ = 0. Thus, 0 = / and is left out 
of our analysis for simplicity. 

In order to understand how such a linearly increas¬ 
ing array of unknown phase shifts may be arranged in 
a practical device, it is useful to consider a specific ex¬ 
ample. Let us suppose that we are to use the QuFTI 
as an optical magnetometer. We consider an interfero¬ 
metric magnetometer of the type discussed in [37] where 
each of the sensing modes of the QuFTI contains a gas 
cell of Rubidium prepared in a state of electromagnet- 
ically induced transparency whereby a photon passing 
through the cell at the point of zero absorption in the 
electromagnetically induced transparency spectrum ac¬ 
quires a phase shift that is proportional to the product 
of an applied uniform (but unknown) magnetic field and 
the length of the cell. We assume that the field is uniform 
across the QuFTI, as would be the case if the entire in¬ 
terferometer was constructed on an all optical chip and 
the field gradient across the chip were negligible. Since we 
are carrying out local phase measurements (not global) 
we are not interested in the magnitude of the magnetic 
field but wish to know if the field changes and if so by how 
much. (Often we are interested in if the field is oscillat¬ 
ing and with what frequency.) Neglecting other sources 
of noise then in an ordinary Mach-Zehnder interferome¬ 
ter this limit would be set by the photon shotnoise limit. 
To construct a QuFTI with the linear cascade of phase 
shifters, as shown in Fig. 1, we simply increase the length 
of the cell by integer amounts in each mode. The first cell 
has length L, the second length 21/, and so forth. This 
will then give us the linearly increasing configuration of 
unknown phase shifts required for the QuFTI to beat the 
SNL. 

One might question why one would employ a phase 
gradient rather than just a single phase. Investigation 
into using a single phase in 4> indicates that this yields no 
benefit. We conjecture that this is because the number 
of paths interrogating a phase in a single mode is not 
superexponential as is the case when a phase gradient is 
employed. 

The interferometer may always be constructed effi¬ 
ciently following the protocol of Reck et al. [38], who 
showed that an n x n linear optics interferometer may 
be constructed from 0{n^) linear optical elements (beam¬ 
splitters and phase-shifters), and the algorithm for deter¬ 
mining the circuit has runtime polynomial in n. Thus, an 
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experimental implementation of our protocol may always 
be efficiently realized. 

The input state to the device is i.e. single pho¬ 

tons inputed in each mode. If ip = 0 then = / and thus 
U = V ' I' = I. In this instance, the output state is 

exactly equal to the input state, Thus, if we de¬ 

fine P as the coincidence probability of measuring one 
photon in each mode at the output, then P = 1 when 
if = 0. When (p ^ 0^ in general P < 1. Thus, intuitively, 
we anticipate that P{p) will act as a witness for p. 

In the protocol, assuming a lossless device, no measure¬ 
ment events are discarded. Upon repeating the protocol 
many times, let x be the number of measurement out¬ 
comes with exactly one photon per mode, and y be the 
number of measurement outcomes without exactly one 
photon per mode. Then P is calculated ns P = x/{x-\-y). 
Thus, all measurement outcomes contribute to the signal 
and none are discarded. Note that, due to preservation 
of photon-number and the fact that we are considering 
the anti-bunched outcome, P{p) may be experimentally 
determined using non-number-resolving detectors if the 
device is lossless. If the device is assumed to be lossy, 
then number-resolving detectors would be necessary to 
distinguish between an error outcome and one in which 
more than one photon exits the same mode. The circuit 
for the architecture is shown in Fig. 1. 


|i> 

|i> 

|i) 


u 



FIG. 1: Architecture of the quantum Fourier transform 
interferometer (QuFTI) for metrology using single-photon 
states. The input state comprises n single photons, 

The state evolves via the passive linear optics unitary 
^ = y . i). 0 . ^Pere y is the quantum Fourier trans¬ 
form, $ is an unknown, linear phase gradient, and © is a 
reference phase gradient used for calibration. At the output 
we perform a coincidence photodetection projecting on ex¬ 
actly one photon per output mode, measuring the observable 
O = (|1)(1|)^’^, which, over many measurements, yields the 
probability distribution P{p) that acts as a witness for the 
unknown phase p. 

The state at the output to the device is a highly path- 
entangled superposition of terms, which grows su- 

perexponentially with n. This corresponds to the number 
of ways to add n non-negative integers whose sum is n, or 
equivalently, the number of ways to put n indistinguish¬ 
able balls into n distinguishable boxes. We conjecture 
that this superexponential path-entanglement yields im¬ 


proved phase-sensitivity as the paths query the phases a 
superexponential number of times. 

The observable being measured is the projection onto 
the state with exactly one photon per output mode, 
O = (|1)(1|)^’^. Thus, (O) = (O^) = P. And, the phase- 
sensitivity estimator reduces to. 


Ap 


dP 

d(f> 


(6) 


Following the result of [39] , P is related to the perma¬ 
nent of P as. 


P=|Per([/)|t (7) 

Here the permanent of the full n x n matrix is computed, 
since exactly one photon is going into and out of every 
mode. This is unlike the boson-sampling protocol [25] 
where permanents of sub-matrices are computed. 

We will now examine the structure of this permanent. 
The matrix form for the n-mode unitary is given by. 


j,k / 2iTT(j-k) \ ’ 

n I e n — j 

as derived in App. A. Taking the permanent of this ma¬ 
trix is challenging as calculating permanents are in gen¬ 
eral #p -hard. However, based on calculating Per(P^’^^) 
for small n, we observe the empirical pattern. 


Per(P^^^) 


1 

^n-l 


n—1 

n 


i=i 


+ n-j , 


( 9 ) 


as conjectured in App. B. This analytic pattern we ob¬ 
serve is not a proof of the permanent, but an empiri¬ 
cal pattern — a conjecture — that has been verified by 
brute force to be correct up to n = 25. Although we don’t 
have a proof beyond that point, n = 25 is well beyond 
what will be experimentally viable in the near future, and 
thus the pattern we observe is sufficient for experimen¬ 
tally enabling super-sensitive metrology with technology 
available in the foreseeable future. 

Following as a corollary to the previous conjecture, the 
coincidence probability of measuring one photon in each 
mode is. 


P = 


Pev{P^^) 

^ n—1 

-,2n—2 n[« „(j)cos(n99)+6„(j)], (10) 


j = l 


as shown in App. C, where 


«n(j) = 2j{n-j), 

bn{j) = p-2jn + 2f. (11) 

The dependence of P on n and p is shown in Fig. 2. 
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FIG. 2: Coincidence photodetection probability P against the 
unknown phase (p and the number of photons and modes n. 
As n increases, the dependence oi P on p increases, resulting 
in improved phase-sensitivity. 



FIG. 3: Phase-sensitivity /\p against the number of photons n 
(red circles). The shotnoise limit of 1/(black squares) and 
Heisenberg limit of 1/iV (orange triangles) are shown for com¬ 
parison. The QuFTI exhibits phase-sensitivity significantly 
better than the shotnoise limit, and only slightly worse than 
the Heisenberg limit. 


It then follows that, 


dP 

dp 


n—1 

GITII'D/'/OH \ 

{j ) 

|sm(^n(/:?;| 2_^ 
i=i 

an{j)cos{n(p) + bn{j) 


( 12 ) 


as shown in App. D. 

Finally, we wish to establish the scaling of /\p. With 
a small p approximation {sm{p) ^ p^ cos{p) ~ 1 — |(/9^) 
we find. 


Ap 



(13) 


as shown in App. E. Thus, the phase sensitivity scales as 
Ap = as shown in Fig. 3. 

We would like to compare the performance of our 
QuFTI to an equivalent multimode interferometer base¬ 
line for which we will construct the shotnoise limit (SNL) 
and Heisenberg limit (HL). This is a subtle compari¬ 
son, due to the linearly increasing unknown phase-shifts, 
{0, (/9,..., (n — l)p}^ that the QuFTI requires to operate. 
The mathematical relation is shown in Fig. 3, where we 
have converted the number of resources, N, to the num¬ 
ber of photons, n. There is disagreement on how such 
resources should be counted. This is the method, which 
we call Ordinal Resource Counting (ORC), that we feel 
most fairly counts our resources. A more detailed sup¬ 
porting discussion can be found in App. F. 

While computing the sensitivity using the standard er¬ 
ror propagation formula of Eq. 1 provides clear evidence 
that our scheme does indeed beat the SNL, it would be in¬ 
structive to carry out a calculation of the quantum Eisher 
information and thereby provide the quantum Cramer- 
Rao bound, which would be a true measure of the best 


performance of this scheme possible, according to the 
laws of quantum theory. However, due to the need to 
compute the permanent of large matrices with complex 
entries, this calculation currently remains intractable. We 
will continue to investigate such a computation for a fu¬ 
ture work. In general, analytic solutions to matrix per¬ 
manents are not possible. In this instance, the analytic 
result is facilitated by the specific structure of the QuETI 
unitary. Other inhomogeneous phase gradients may yield 
analytic results, but we leave this for future work. 

In App. G we discuss the efficiency of the QuETI pro¬ 
tocol and in App. H we analyse dephasing, which is a 
source of decoherence, and find that the QuETI proto¬ 
col is far more robust against dephasing than the NOON 
state is. 

We have shown that a passive linear optics network fed 
with single-photon Eock states may implement quantum 
metrology with phase-sensitivity that beats the shotnoise 
limit. Unlike other schemes that employ exotic states 
such as NOON states, which are notoriously difficult to 
prepare, single-photon states may be readily prepared in 
the laboratory using present-day technology. This new 
approach to metrology via easy-to-prepare single-photon 
states and disjoint photodetection provides a road to¬ 
wards improved quantum metrology with frugal physical 
resources. 
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One can see the pattern that emerges is of the form: 


Per(C/(")) 


1 


n—1 

n 


i=i 


+ n-j , 


(Bl) 


Appendix A: Proof of 


Beginning from Eq. Al and setting 0 = /, 


which is Eq. 9 that we set out to show. This equation has 
been verified analytically up to n = 16 and up to n = 25 
numerically.. 
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Lm=l 


^ ^—2ijl7r/n ^2imk7rln 
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n ^ 

1=1 




-1 

1 2il(k-j)7l 


= -E 
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-\-i{l — l)ip 


= e 
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2i(k — j) 7^ I V-> 2i(k — 7)7r , • 7 

- y^{e ri ^ 

^ z=o 


From the geometric series, it follows. 
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Appendix C: Calculation of P 

Assuming our conjecture in Eq. 9 holds, we can com¬ 
pute the coincidence probability of measuring one photon 
in each mode at the output. 


P = |Perm(f7(”))|^ 

^ n—1 

;^n "-j) 

i=i 

^ n—1 

i=i 

^ n—1 

',2n—2 n| jcos{nip) + ijsm{n(p) P n — j 

i=i 

^ n—1 

In—2 n| jcos{nip) {n — j) pi jsm{nip) 


-,2n- 


i=i 


Re 


Im 


(Cl) 


Invoking the property that = Re(^)^ + Im( 2 :)^, where 
2: G C, 


which is what we set out to prove, which is Eq. 8 that 
we set out to prove, where the last line follows from the 
geometric series. 


Appendix B: Conjecture for the Analytic Form of 
Per(E(^)) 


Our goal is to find the analytic form for Per(L'^’^^) 
where is as in Eq. Al. We can perform a brute force 
calculation to obtain the analytic form for small n. Doing 
so up to n = 6 yields: 


n 

PerCt/f"-)) 

1 

1 

2 

cos(iji) 

3 

1 (2 + (1 + 26®’“^) 

4 

^ (1 + (3 + (1 + 

5 

^ (4 + (3 + 26®"''') (2 + 36®*“^) (l + 

6 

^ (1 + (2 + e®"''') (5 + e®**^) (l + 26®'“^) (l + Se®"''') 


It — ± 

^ = 2n-2 n \{jcos{nip) + {n - jp + fsin^inip) 

n , i 

j=i 

^ n —1 

= ^^2^132 n [ yCOS^ {nip) + {rvp) 

i=i ' 2^2 ' 

+ ‘^3{n - j)cos{rvp) + (n - jf 

= n [-1^ + “ j)cos{rup) + (n - j)^ 

n , i 

j=i 

^ n —1 

= ^2n-2 n i~ cos(n<^) + y - 2 jn + 2 f^ 

a„(i) 6„(i) 

^ n—1 

= n [anij)cos{nip) + b„{j) , 

^ i=i 

(C2) 


which is Eq. 10 that we set out to show. 
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Appendix D: Calculation of ||^| 

From Eq. C2, exploiting the logarithm product rule. 

ln(P) = ln(^-^^+ln | [« n{j)cOs{n(f) +bn{j) 

s_✓ \i=i 

c 

n—1 

= c + E ln[an(j)cos(n(y?) + bn{j) , (Dl) 

i=i 

where C is a constant. Now the derivative becomes, 

1 dP nan{j)sm{n(p) 

P dip ^ an{j)cos{nip) + bn{j) 


^ ^ -nPsin(nM)V__ 


Thus, 

dP 


dip 


n—1 

= nP\sm{nip)\ ^ 
i=i 


{j ) 


ttn(i)cos(n(/9) + 6n(j) 


(D2) 


, (D3) 


where /c(n) = in(n — l)(n + 1) > 0 V n > 1. From Eq. 


compute and to be. 



' Si 1 — 2k{n)p 

(E3) 

dP 

dp 

= 2k{n)\<f\, 

(E4) 


where we have again dropped terms above order ip‘^. Us¬ 
ing Eq. 6 the phase sensitivity Aip in the small angle 
regime is. 


Aip = 




dP 

d(f> 


— k{n)ip‘^^ ~ ~ 2k{n)ip‘^^ 


2k{n)\ip\ 


VHri)ip‘^ 

2k{n)\ip\ 

1 

2^/Hr^ 


y 2(n — l)n{n + 1) ’ 

which is Eq. 13 that we set out to show. 


(E5) 


which is Eq. 12 that we set out to show. 


Appendix E: Calculation of Ap in the small angle 
approximat ion 


We wish to compute Ap in the limit that np <C 1. 
Then P in the small angle regime of Eq. 10 becomes. 


1 r 1 


i=i 

n—1 


n—1 


^2n—2 J.J. 


i=i 


- {nj Pp)n^p‘^ 


= n 


1=1 


(El) 


where cos(n(^) is expanded to the first nonconstant term 
in its Taylor series. This product has the form of a bi¬ 
nomial expansion. Dropping terms above order p‘^ ^ P re¬ 
duces to. 


n—1 

1=1 

= 1 — p‘^ -(n —l)n(n + l) 
= 1 — k{n)p‘^, 


(E2) 


Appendix F: Discussion of Ordinal Resource 
Counting (ORC) 

We would like to compare the performance of our 
QuFTI to an equivalent multimode interferometer base¬ 
line for which we will construct the shotnoise limit (SNL) 
and Heisenberg limit (HL). This is a subtle compari¬ 
son, due to the linearly increasing unknown phase-shifts, 
{0, (^,..., (n — l)ip}, that the QuFTI requires to operate. 
There is a long and muddled history of increasing the in¬ 
terrogation time (or here length) of the probe particles 
with the unknown phase-shift followed by an incorrect 
reckoning of the true resources. Here we shall use a pro¬ 
tocol we call Ordinal Resource Counting (ORC) whereby 
all resources, such as number of ‘calls’ to the phase-shifter 
p, are converted to the ‘currency’ of the resource that is 
most precious to us, namely photon-number. We do this 
as follows. 

First we must construct a multimode interferome¬ 
ter with n photon inputs that provides the baseline if 
the photons remain uncorrelated and the number-path 
entanglement remains minimal. Such a comparator is 
shown in Fig. 4, and consists of n, two-mode Mach- 
Zehnder Interferometers (MZI) in a vertical cascade, fed 
with single-photon inputs, with the same linearly increas¬ 
ing unknown phase-shift sequence as the QuFTI. Since 
the MZIs are disconnected, the number-path entangle¬ 
ment remains constant and minimal, and of the form 
(|1, 0) -t- |0, l))/\/2 inside each MZI. 
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ters, with a linearly increasing phase gradient. This system 
has the same configuration of phases as the QuFTI, but the 
photons are not allowed to interfere, and thus has minimal 
number-path entanglement. 


|i)« 



FIG. 5: Noting that a single photon interrogating a phase- 
shift of n(f is equivalent to n independent interrogations of 
9?, Fig. 4 can be represented in terms of the resource of pho¬ 
tons as shown here. Here means that j independent (i.e 
distinguishable) photons have been prepared. 


Now to convert the linearly increasing interrogation 
lengths of the unknown phase-shifts, we note that a sin¬ 
gle photon interrogating a phase-shift of say 2(p is equiv¬ 
alent to a single photon interrogating a single phase-shift 
(f twice, which is in turn equivalent to two uncorrelated 
photons entering the same port of the MZI containing 
a single phase-shift of cp. In this way we may convert 
‘number of interrogations of the phase-shifter’ into the 
currency of ‘number of photons’ to carry out a fair reck¬ 
oning of the resources. Following this logic we are led 
to Fig. 5 showing a cascade of MZIs where the linearly 
increasing phase-shifters are replaced with a single phase- 
shifter of ip and the single photons at the MZI inputs are 
replaced with a linearly increasing number of photons. 
Then the ‘number of interrogations of the phase-shifter’ 
becomes n{n — l)/2, but there is an additional photon 
that is part of the QuFTI resources so our total number 
of resources becomes, 


N 


= 1 + ~ 1 ) 


(Fl) 


Next we note that this cascade of n MZIs in Fig. 5 may 
be replaced with a single MZI, shown in Fig. 6, where 


the input is now an ordinal grouped ranking of the un¬ 
correlated photons following the same pattern as in Fig. 
5. Hence in the configuration in Fig. 6 we have a sin¬ 
gle MZI with vacuum entering the lower port, a stream 
of N uncorrelated photons entering the upper port, and 
a single phase-shifter ip between the beamsplitters. It is 
well-known that for this configuration the sensitivity of 
this system scales as the SNL [3, 40], namely. 


A^snl = ^ = 


1 + 


n(n—1) 


(F2) 


This then provides us a fair reckoning of the SNL to be 
used gauging the performance of the QuFTI. 

Finally, if instead we were to maximally path-number 
entangle these resources into a NOON state of the form 
(|A^, 0) + |0,A^))/v^ (just to the right of the first beam 
splitter but before the phase-shifter) the sensitivity then 
becomes Heisenberg limited. 


Av^hl = ^ = 


1 


1 + 


n(n—1) 


(F3) 


which is a sensitivity known to saturate the Quantum 
Cramer-Rao Bound (CRB) for sensitivity in local phase 

















































FIG. 6: Grouping all the independent interferometers in Fig. 
5 together and including the extra photon from the QuFTI 
model, we obtain a single MZI with 1 + n(n — l)/2 indepen¬ 
dent photons as input. This configuration achieves the shot- 
noise limit, and thus provides a benchmark for comparing our 
QuFTI protocol against the shot noise and Heisenberg limits, 
with photons as the resource being counted. 



n 

FIG. 7: Dephasing for (f = 0.01. The shaded region represents 
the phase sensitivity for the QuFTI where 0 < y < 0.01. 


estimation with N photons [12, 13]. As the CRB is the 
best one may do, according to the laws of quantum me¬ 
chanics, then in this case the HL is optimal. As discussed, 
the performance of the QuFTI falls between the SNL and 
the HL, but with the feature of not having to do anything 
resource intensive such as preparing a high-NOON state. 

Thus the SNL and the HL, computed via this Ordinal 
Resource Counting method, provides the fairest compar¬ 
ison of sensitivity performance of the QuFTI with such 
ambiguities such as how to handle ‘number of calls to the 
phase-shifter’ removed by replacing such a notion with 
‘number of photons’ inputted into the interferometer. 


Appendix G: Efficiency 

In the presence of inefficient photon sources and photo¬ 
detectors the success probability of the protocol will drop 
exponentially with the number of photons. Specifically, 
if r]s and r]d are the source and detection efficiencies 
respectively, the success probability of the protocol is 
V = iVsVd)^• Current cutting edge transition edge detec¬ 
tors operate at 98% efficiency, with negligible dark count 
[41]. SPDC sources are the standard photon-source tech¬ 
nology but they are non-deterministic. However, there 
are techniques that can greatly improve the heralding ef¬ 
ficiency up to 42% at 2.1 MHz [42]. Also, other source 
technologies, such as quantum dot sources are becoming 
viable with efficiencies also up to 42% [43]. For n = 10, 
which is already well beyond current experiments, this 
yields r] = (0.98 * 0.42)^^ ^ 0.00014, which is about 300 
successful experimental runs per second when operating 
with 2.1 MHz sources. 


Appendix H: Dephasing 

A form of decoherence to consider is dephasing. De¬ 
phasing in our work may be modelled with the result of 
Bardhan et al. [44], whereby dephasing occurs on each 
mode separately. When considering our example of a 


magnetometer, dephasing would occur in the magnetic 
field cells where atomic fluctuations may occur that differ 
between cells. In the rest of the interferometer, dephas¬ 
ing can be made very close to zero, particularly on an all 
optical chip. 

To model dephasing we investigate a random phase 
shift Ax added to each mode separately. Ax is a Gaussian 
random variable of zero mean but nonzero second order 
moment. The phase shift in the jth mode then becomes, 

_ ^±ijcp^±ijAx 

= 6='=*-^'^ (^1 ± ijAx - ± ... y (HI) 

Using (Ax) = 0, (Ax^) % 0, and that Ax ^ 0 we sim¬ 
plify this to be, 

« (H2) 

The signal P in Eq. 10 from our work then changes in 
the presence of dephasing. The dependence that P has 



FIG. 8: The effect of dephasing on the NOON state and 
QuFTI where (p = 0.01, x = 0.005. The NOON state is plot¬ 
ted with respect to N for fair resource counting. 
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on the unknown phase Lp does not depend on the mode 
number j. Then the term that depends on ip becomes, 

cos(n</.) = i 

= cos(n0)e“^’^ (H3) 


The factor e can be absorbed into an{j) so that 

the derivation of |^| in Eq. 12 is identical. Using this 
result we numerically plot the phase sensitivity with de¬ 
phasing in Fig. 7. 


Using this substitution P becomes, 
|2 


P = 


Per(U^^^) 


^ n—1 

= n [«n(i)cos(n?!.)e"5""'^^" +&n(j) ;H4) 


j=l 


In order to meaningfully analyze the dephased sensi¬ 
tivity, we would like to compare with other well known 
metrological schemes. In Fig. 8, we compare the QuFTI 
to the NOON state (with N input photons for a fair re¬ 
source comparison) and see that the QuFTI is far more 
robust against dephasing. 
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